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ABSTRACT 


A  photoreaction-diffusion  model  for  the  solar  degradation  of  thermal 
control  materials  is  mathematically  formulated  and  analyzed.  The 
governing  differential  equation  for  the  concentration  of  defect  specie 
created  in  the  material  by  radiation  is  derived,  and  an  analytical  solu¬ 
tion  is  obtained  for  appropriate  boundary  conditions.  From  this  solution 
for  the  concentration  Ni  of  the  ith-type  defects,  an  analytical  expres¬ 
sion  is  derived  for  the  number  Qj  of  the  ith-type  defects  present  in  a 
material  layer  of  depth  x  and  unit  irradiated  area.  These  analytical 
relations  are  used  to  predict  the  radiation-induced  increase  in  the  ab¬ 
sorption  coefficient  and  absorption  optical  thickness  of  a  hypothetical 
thermal  control  material  exposed  to  damaging  ultraviolet  irradiance 
equal  to  that  of  the  sun.  The  solutions  for  Ni  and  Qi  are  also  used  to 
better  define  the  environmental  parameters  which  must  be  controlled 
during  solar  radiation  testing  of  thermal  control  coatings.  For  example, 
the  time -irradiance  reciprocity  principle  often  employed  in  accelerated 
solar  testing  is  investigated  by  use  of  the  aforementioned  solutions  and 
is  found  to  be  invalid  for  materials  which  degrade  according  to  the 
photoreaction -diffusion  model. 
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NOMENCLATURE 

D  Diffusion  coefficient 

erfc  Complementary  error  function 

F_j_  Function  defined  in  Eq.  (8) 

G_i  Function  defined  in  Eq.  (13) 

I  Radiant  energy  flux  in  photons /sec-cm^ 

k  Quantum  yield,  net  number  of  defects  created  per  absorbed 

photon 

N  Defect  concentration 

n  Number  of  different  types  of  defect  specie 

P  Pressure 

Q  Number  of  given  type  of  defects  present  in  a  layer  of  depth  x 

and  unit  irradiated  area 

R  Reflectance 

T  Temperature 

t  Irradiation  time 

x  Local  depth  in  material  (see  Fig.  1) 

a  Absorptance 

7  Wavelength  of  damaging  radiation 

A  Change  in  an  associated  quantity 

e  Absorption  coefficient 

X  Wavelength  of  radiation 

aa  Absorption  cross  section 

t  Absorption  optical  thickness 
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SECTION  I 
INTRODUCTION 


In  several  previous  studies  of  solar  ultraviolet  radiation  effects  on 
the  zinc  oxide  (ZnO)  pigments  of  thermal  control  coatings  (Refs.  1 
through  3),  the  investigators  postulated  that  a  photoreaction-diffusion 
damage  mechanism  causes  the  changes  induced  in  the  optical  properties 
of  the  zinc  oxide.  They  also  suggested  that  the  same  type  of  mechanism 
causes  the  ultraviolet  (UV) -induced  damage  observed  for  similar  thermal 
control  materials  such  as  titanium  dioxide  (Ti02)  and  tin  oxide  (SnC>2). 

In  the  work  presented  here,  a  photoreaction-diffusion  model  for  the  solar 
degradation  of  thermal  control  materials  of  this  type  is  mathematically 
formulated  and  analyzed.  That  was  the  main  objective  of  this  study. 
There  were  several  secondary  objectives.  One  was  to  deduce  a  mathe¬ 
matical  model  which,  when  completely  formulated,  might  permit  the 
theoretical  prediction  of  solar  absorptance  change  of  spacecraft  thermal 
control  surfaces  as  a  function  of  irradiation  time,  irradiance  level,  and 
surface  temperature.  An  example  application  of  this  model  is  given 
and  used  to  predict  the  induced  absorption  coefficient  and  the  increase 
in  absorption  optical  thickness  (Ref.  4)  of  a  hypothetical  thermal  con¬ 
trol  material  exposed  to  damaging  UV  irradiance  equal  to  that  of  the 
sun.  Another  objective  was  to  mathematically  define  the  space  environ¬ 
ment  conditions  which  should  be  simulated  in  laboratory  tests  of  solar 
degradation  effects  on  thermal  control  coatings.  Among  the  environ¬ 
mental  effects  considered  are  the  magnitude  and  spectral  distribution 
of  the  damaging  irradiance.  For  example,  a  theoretical  investigation 
was  made  of  the  time -irradiance  reciprocity  principle  used  in  acceler¬ 
ated  solar  radiation  testing. 


SECTION  II 

FORMULATION  OF  MATHEMATICAL  MODEL  FOR  PHOTOREACTION- 
DIFFUSION  DAMAGE  MECHANISM 


2.1  PHYSICAL  MODEL  OF  DAMAGE  MECHANISM 

The  photoreaction-diffusion  damage  mechanism  postulated  for  the 
solar  degradation  of  thermal  control  materials  basically  assumes  that 
ultraviolet  photons  absorbed  near  the  surface  of  the  material  possess 
sufficient  energy  to  produce  defects  by  photochemical  reaction.  The 
defects  are  created  within  the  penetration  depth  of  the  damaging  irradia¬ 
tion  and  then  diffuse  into  the  material  bulk  as  the  generation  process 
continues  at  a  decreasing  rate  near  the  surface.  These  defects  will 
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absorb  radiation  in  a  spectrally  dependent  manner  and  hence  may  cause 
a  significant  increase  in  the  parent  material  absorption  properties  at 
wavelengths  in  the  visible  and  infrared.  It  is  assumed  that  the  change  in 
monochromatic  absorptance  at  these  wavelengths  is  a  function  of  the 
defect  concentration  in  the  parent  material. 


2.2  MATHEMATICAL  MODEL  AND  DERIVATION  OF  GOVERNING  EQUATIONS 

A  mathematical  model  for  the  UV  degradation  of  thermal  control 
materials  may  be  formulated  in  the  following  manner.  The  material  is 
assumed  to  be  in  a  low  pressure  vacuum  environment.  Its  surface,  as 
shown  in  Fig.  1  (Appendix),  is  exposed  to  damaging  irradiation  of  wave¬ 
lengths  lying  within  the  differential  wavelength  intervals  7. ~  The  spectral 
irradiance  at  wavelength  7  lying  within  dy  is  denoted  by  IQj  7  and  the  total 
ir radiance  with  d7  is  given  by  IQj  -ydy.  Most  of  this  incident  radiant 
energy  flux  will  be  absorbed  by  the  material  within  a  depth  of  less  than 
0.  1  micron  from  the  surface  and  there  generate  by  photoreaction 
ith-type  defects  per  absorbed  photon.  Figure  1  shows  a  cross  section 
of  this  region  of  the  material  with  a  differential  element  dx  located  at  a 
distance  x  from  the  surface  (x  =  0).  The  concentration  of  the  ith  defect 
specie  produced  by  damaging  irradiation  of  wavelength  7  and  now  present 
in  dx  at  time  t  is  denoted  by  ydy.  Hence,  the  symbol  T  can  be  in¬ 
terpreted  as  the  concentration  of  the  ith-type  defects  created  per  wave¬ 
length  interval  dy  by  damaging  radiation  of  wavelength  7.  The  net  radiant 
energy  flux  downward  at  depth  x  is  shown  in  Fig.  1  and  denoted  by  Ix  ydy. 

9Ni 

Figure  1  also  shows  the  diffusion  fluxes  -Di  ~ —  dy  which  exist  in  the 

material  because  of  the  presence  of  specie  i.  The  symbol  Dj  denotes  the 
diffusion  coefficient  of  the  ith  defect  specie  in  the  parent  material.  It  is 
assumed  that  there  are  no  interactions  between  defects. 

To  derive  the  governing  differential  equations  for  the  defect  species 
concentrations,  the  law  of  mass  conservation  is  assumed  for  the  ith 
defect  specie  present  in  the  differential  element  dx  shown  in  Fig.  1.  This 
yields  the  following  partial  differential  equation 


where  kiy  is  designated  as  the  ith  defect  quantum  yield  (Ref.  5)  of  the 
material  for  damaging  radiation  of  wavelength  7.  From  this  equation  it 
can  be  seen  that  an  analytical  expression  for  Ix>  y  is  needed.  For  a 
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collimated  radiant  flux  incident  normally  on  the  surface  as  shown  in 
Fig.  1,*  the  analytical  relation  for  IXj  y,  assuming  internal  scattering 
to  be  negligible  at  wavelength  y,  is 


Hence, 


Ix.y  =  (1  -  Ry)I0,yexP  (“  «y*) 


dlx,y 

dx 


-  (1  -  Ry)  fy  I o,y  exp  (-fyX) 


(2) 

(3) 


where  Ry  and  €y  are,  respectively,  the  reflectance  and  linear  absorp¬ 
tion  coefficient  for  the  parent  material  at  the  wavelength  of  the  damaging 
radiation.  Thus,  from  Eq.  (1)  with  Di  independent  of  x  and  Eq.  (3),  the 
governing  differential  equation  for  the  concentration  of  defect  specie  i  is 


dNi.y 

di 


kjy(l  -  Ry)  tv  I0,y  exp  (- Cyx) 


D. 


d2  Ni,y 

dx2 


(4) 


If  n  different  types  of  noninteracting  defect  specie  are  generated  by 
damaging  radiation  of  wavelength  y,  there  will  be  n  different  equations 
of  the  form  of  Eq.  (4). 


2.3  FACTORS  POTENTIALLY  INFLUENCING  PHYSICAL 
PARAMETERS  l0#  y,  Dj.  (y,  R>/(  AND  kiy 

As  noted  by  the  subscript  y,  the  physical  parameters  I0,  y,  £y,  R-y, 
and  appearing  in  Eq.  (4)  can  all  be  a  function  of  the  wavelength  of 
the  damaging  radiation.  These  parameters  along  with  Dj  may  also  be  a 
function  of  other  factors.  For  example,  the  spectral  irradiance  I0,  y 
might  be  a  function  of  time,  t.  Potentially  influencing  factors  for  the 
other  parameters  are  described  below. 

First,  consider  Di,  the  diffusion  coefficient  of  the  ith  defect  specie. 
In  general,  Di  can  depend  on  the  temperature  of  the  parent  material, 
the  concentration  of  the  ith-type  specie,  and  also  the  concentrations  of 
any  other  types  of  specie  which  may  be  present.  However,  the  depend¬ 
ence  of  Di  on  concentration  is  usually  negligible  except  at  very  high  con¬ 
centration  levels  (Ref.  6). 

Next,  consider  the  absorption  coefficient  and  reflectance  of  the 
parent  material  at  the  wavelength  of  the  damaging  radiation,  €y  and  Ry, 
respectively.  In  general,  it  is  felt  that  €y  will  not  vary  greatly  with 


*In  Fig.  1,  it  is  assumed  that  the  penetration  depth  of  the  damaging 
irradiation  is  much  less  than  the  material  thickness. 
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location  in  the  thermal  control  material  and  hence  can  be  considered  a 
constant  with  respect  to  x.  This  assumption  is  employed  throughout  the 
analysis.  In  addition,  it  is  assumed  that  e-y  is  not  a  function  of  time,  t. 
This  assumption  physically  means  that  although  defects  generated  by 
damaging  irradiation  of  wavelength  7  can  absorb  radiation  at  that  same 
wavelength,  the  resulting  increase  in  absorption  at  7  is  negligible  rela¬ 
tive  to  the  initial  absorption  at  7.  An  illustrative  problem  presented 
later  in  the  text  will  prove  that  this  assumption  is  definitely  valid.  As 
for  R-y,  the  reflectance  of  the  material  at  the  wavelength  of  the  damaging 
radiation,  it  also  is  assumed  constant  with  respect  to  time,  t. 

Finally,  it  is  necessary  to  consider  kj/yj  the  ith  defect  quantum 
yield.  As  noted  earlier,  ki-y  is  a  function  of  the  wavelength  of  the 
damaging  radiation  (Ref.  5).  It  also  can  be  a  function  of  the  radiant 
exposure  of  the  surface.  For  thermal  control  materials  of  the  type 
mentioned  previously,  ki-y  is  likely  to  decrease  rapidly  as  the  radiant 
exposure  of  the  surface  increases  (Ref.  7).  It  is  speculated  that  kj-y 
is  also  a  function  of  the  irradiance  level  of  the  damaging  radiation 
(Ref.  8).  In  addition,  the  temperature  of  the  material  can  affect  k^-y 
through  its  influence  on  secondary  chemical  processes  taking  place  in 
the  medium  (Ref.  8).  It  is  also  possible  that  the  pressure  of  the  vacuum 
environment  in  which  the  material  is  located  could  affect  ki-y.  This 
would  happen  if  the  pressure  were  not  maintained  at  a  low  enough  level 
to  ensure  that  the  concentration  of  any  gases  present  was  insufficient  to 
influence  the  secondary  chemical  processes  occurring  in  the  material 
and  on  its  surface  (Ref.  7). 


SECTION  III 

AN  ANALYTICAL  EXPRESSION  FOR  THE  CONCENTRATION  OF 
ith-TYPE  DEFECT  SPECIE 


An  analytical  solution  has  been  derived  to  Eq.  (4)  for  I0j  j  constant 
with  time  and  D*  independent  of  concentration.  This  solution  was  ob¬ 
tained  under  the  following  assumptions:  (1)  Initially  there  are  no 
radiation-induced  defects  present  in  the  material,  (2)  once  the  irradia¬ 
tion  begins,  none  of  the  defects  created  by  it  will  be  lost  through  the 
irradiated  surface  (x  =  0),  and  (3)  none  of  the  radiation-induced  defects 
will  penetrate  to  the  material  boundary  opposite  the  irradiated  surface. 
Mathematically,  these  assumptions  may  be  expressed  as  the  following 
initial  and  boundary  conditions  on  Eq.  (4): 
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t  =  0,  Nj.y  =  0  for  x  >  0 


( 


t  >  0,  <  (9x 


s  ax 

W  - 


=  0 

X  =  0 

0  as  x  -*  °= 


(5) 


The  solution  to  Eq.  (4),  subject  to  the  conditions  stated  in  Eq.  (5),  has 
been  obtained  by  operational  techniques  for  kiy  being  the  following 
function*  of  radiant  exposure, 


kjy  -  kiyr 


/4y_\ 

Vr,ytr/ 


(6) 


This  solution  to  Eqs.  (4),  (5),  and  (6)  takes  the  form 


N‘.r  ■  (-f)  0  w  V  F-  [2‘y •  J7^r]  (7) 


where,  noting  eyx  =  2fy(Djt)!4  . 

r  * 

E-,  2fy  (Djt) 


2(Djt) 


2(Djt)K 
=  erfc 


(8) 


-exp(«y2'  D £t  -  tyx),  .  erfc  [fy  (Djt)^l 
+  •  exp  (fy1  Djt  -  fyx).  .  erfc  [ey  (Djt)  ^  -  - - — ^-1 

2  l*  2(D;t)^  J 

-  ^  •  exp  (fy*  Djt  +  fyx)  .  erfc  [fy(Djt)^  +  — - — u-1 

L  2(Dit)/4  J 


Figure  2  presents  this  solution  for  y  in  the  dimensionless  form  F_j 

with  2^.  ^1/2  ^a^en  as  t^le  independent  variable  and  26y(Djt)i/  as  a 

parameter.  Note  that  F_i  is  essentially  the  same  for  all  values  of 
2e7(D1t)1/2  >  10. 

An  analysis  of  Eq.  (7)  indicates  that  it  can  be  used  to  predict  analyt¬ 
ically  the  x -distribution  of  y  as  a  function  of  time,  t,  and  the  spectral 
irradiance,  IGj  y.  Thus,  yiy,  the  concentration  of  the  ith  defect  specie 
generated  by  damaging  radiation  lying  within  the  small  wavelength  interval 
dy,  can  be  predicted  as  a  function  of  the  same  variables.  For  the  case 


*The  subscript  r  on  the  symbols  in  Eq.  (6)  denotes  reference  values 
of  the  associated  physical  quantities. 
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where  the  material  is  irradiated  with  total  radiation,  the  concentration 
of  the  ith  defect  specie,  Ni,  must  be  obtained  by  integrating  Ni>  y  of 
Eq.  (7)  over  all  damaging  wavelengths  since  photons  at  other  wave¬ 
lengths  besides  7  may  also  produce  defect  specie  of  type  i.  The  result¬ 
ing  expression  for  Ni  is 


yicn 
(l  - 


Ry) 


kiyr  dr,y  Io,y)  F— i  |  2«y(9il) 


4* 


2  (D;l) 


3 


dy 


(9) 


where  7^  is  the  largest  wavelength  at  which  incident  photons  will  produce 
defect  specie  of  type  i.  It  is  indicated  that  Yj  can  be  a  function  of  the 
material  temperature,  T.  This  is  because  the  wavelength  at  which  the 
material  absorption  edge  occurs  may  shift  with  temperature  (Ref.  3). 
Note  that  the  concentration  of  defects  at  the  surface  can  be  obtained 
directly  from  Eq.  (9)  and  is 


.  *  ~Yi 

Nio  =  J  (1  -  Ry)  kiyr(Ir,y  Io.y)*  F_,  [2fy  (Dit)54,  0]  dy 

1  '  o 

F_,  [2fy  (DjO^.O]  «  1  -  exp  Uy1  Djt)  erfc  [fy  (Djt)5*] 


(10) 


An  analysis  of  Eq.  (9)  indicates  that  it  can  be  used  to  predict  the 
concentration  distribution  of  ith-type  defects  created  in  the  irradiated 
material  by  radiation  of  wavelengths  lying  within  the  wavelength  interval 
(0,  Yj).  This  can  be  done  as  a  function  of  irradiation  time  and  the  magni¬ 
tude  and  spectral  distribution  of  the  damaging  irradiance.  Such  an  ex¬ 
pression  is  useful  because  the  concentration  Ni  determines  the  increase 
in  the  spectral  absorption  coefficient  of  the  material  caused  by  the 
presence  of  ith-type  defect  specie. 


SECTION  IV 

NUMBER  OF  ith-TYPE  DEFECTS  PRESENT  IN  A  LAYER  OF 
DEPTH  x  AND  UNITIRRADIATED  AREA 


The  number  Qi  of  ith-type  defects  present  in  a  layer  of  depth  x  and 
unit  irradiated  area  is  mathematically  related  to  the  specie  concentration 
Ni  by  the  expression 


(11) 
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Thus,  Qi  is  the  integral  of  the  specie  concentration  over  depth  x  into 
the  material.  This  integral  has  been  evaluated  for  Ni  given  by  Eq.  (9) 
and  the  resulting  expression  is 


(1  -  H  y)  kiyc  (I  t,y  1 0,y)A  [2,y  <D,t>\ 


dy  (12) 


where 


-^^77 j  =  exp  (ty2  Dit  -  fyx)  •  i  erfc  [ey(Djt)^] 

+  ('  -  cxp(w)) +  •»* er,c  [777] 


1 

2 


1 

2 


exp  (fyJ  Djt  -  £yx)  •  erfc  j«y  (D;t)^  — 

r  % 

exp  ( ty 1  Dit  +  fyx)  ■  erfc  lfy( D ;t)  + 


(13) 


The  function  G_i  is  displayed  in  Fig.  3  with  l/~2  as  the  independ- 

1/2 

ent  variable  and  2e,y(Djt)  '  taken  as  a  parameter.  Note  that  G_i  can 
also  be  interpreted  as 

q.  j{lrlr/ ni)K  0  ~  M  (lp.y  (14) 


Qi.y 


dx  / 


Equation  (12)  for  Qi  can  be  used  to  predict  the  number  of  ith  defects 
present  in  a  material  layer  of  depth  x  and  unit  irradiated  area  as  a  func¬ 
tion  of  irradiation  time  and  the  magnitude  and  spectral  distribution  of 
the  damaging  irradiance.  Such  an  expression  is  particularly  useful 
since  Qi  strongly  influences  the  increase  in  absorption  optical  thickness 
resulting  from  the  presence  of  ith-type  defects.  This  increase  in  ab¬ 
sorption  optical  thickness  is  one  of  the  most  important  parameters 
governing  the  induced  spectral  absorptance  of  thermal  control  materials. 
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SECTION  V 

RELATION  OF  N;  AND  Qf  TO  THEORETICAL  PREDICTIONS 
OF  SOLAR  ABSORPTANCE  CHANGE 


5.1  GENERAL  DISCUSSION  OF  MATHEMATICAL  MODEL 

In  a  mathematical  model  for  predicting  the  solar  absorptance  change 
of  spacecraft  thermal  control  surfaces  as  a  function  of  irradiation  time 
and  irradiance  level,  it  is  assumed  that  defects  are  generated  in  thermal 
control  materials  by  absorption  of  high  energy  photons.  If  the  defect 
concentration,  Ni,  and  the  number,  Qj,  of  ith-type  defects  per  unit 
irradiated  area  could  be  predicted,  then  knowledge  of  the  basic  absorp¬ 
tion  properties  of  the  defect  would  allow  predictions  of  the  increase  in 
the  material  spectral  absorption  coefficient  and  absorption  optical  thick¬ 
ness.  With  the  increases  in  absorption  coefficient  and  absorption  optical 
thickness  determined,  one  would  then  solve  the  radiative  transport  equa¬ 
tion  with  these  variables  present.  The  resulting  solutions  for  the  in¬ 
tensity  and  radiant  flux  could  be  used  to  predict  the  change  in  the  material 
spectral  absorptance  as  a  function  of  irradiation  time  and  irradiance  level. 

Obtaining  analytical  relations  between  the  irradiation  time  and 
and  Qi  would  be  an  important  step  in  any  model  used  for  theoretically 
predicting  the  change  in  solar  absorptance  of  thermal  control  materials. 
Another  important  step  would  be  to  develop  an  analytical  relation  between 
Ni,  Qj  and  the  spectral  absorptance  increase  of  these  materials. 

5.2  ASSUMED  RELATIONS  BETWEEN  N-,  Q.  AND  THE  INDUCED  ABSORPTION 
PROPERTIES  OF  THERMAL  CONTROL  MATERIALS 

Any  attempt  to  develop  analytical  relations  between  Nj,  Qi  and  the 
increase  in  spectral  absorptance  of  thermal  control  materials  will  first 
require  that  theoretical  expressions  be  established  between  the  induced 
absorption  properties  of  the  material  and  Ni  and  Qi.  This  study  assumes 
that  a  linear  relation  exists  between  the  induced  absorption  coefficient 
and  the  defect  specie  concentration.  Thus,  the  increase  in  the  mono¬ 
chromatic  absorption  coefficient  of  the  parent  material  caused  by  the 
presence  of  ith-type  defect  specie,  Ae^,  is  the  product  of  Nj  and  the  ab¬ 
sorption  cross  section  for  the  ith  defect  specie  at  wavelength  X,  aaix. 

A(\i  =  ffaiXNj  (15) 

where  Ae^i  will  be  referred  to,  for  conciseness,  as  the  induced  mono¬ 
chromatic  absorption  coefficient  of  the  material  for  type  i  defect  specie. 
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When  n  different  types  of  defect  specie  are  present,  the  total  induced 
absorption  coefficient  of  the  material,  at  wavelength  X,  Ae^,  is  assumed 
to  be  the  linear  sum  of  the  effect  of  each  type  of  specie,  or 

n  n 

A'A  =  2]  A‘Ai  =  T.  aaiA  Nj  (16) 

i=i  i=i 

This  equation  and  the  assumption  on  which  it  is  based  ignore  any  inter¬ 
action  effects  that  might  occur  between  different  defect  species. 


Once  a  relation  between  Ae^  and  is  specified,  it  is  possible  to 
establish  a  connection  between  the  increase  in  monochromatic  absorp¬ 
tion  optical  thickness,  At^,  and  N^.  For  conciseness,  Ar^  will  be 
referred  to  as  the  induced  absorption  optical  thickness  of  the  material 
at  wavelength  X.  An  analytical  expression  between  Ar^  and  Nj  is  ob¬ 
tained  from  Eq.  (16)  and  the  definition  of  At^, 

Jx  /»* 

A<rAd*  =  £  0-aiXj  Nidx'  =  £  ^aixQi  (17) 

o  j  =  l  0  i  =  1 


Equation  (17)  illustrates  the  relation  between  At^  and  Ni,  or  Qi.  This 
relation,  of  course,  implicitly  contains  the  assumptions  on  which 
Eqs.  (15)  and  (16)  were  based.  Note  that  the  induced  absorption  optical 
thickness  resulting  from  the  presence  of  an  ith-type  defect  specie,  At^j, 
is  given  by  the  expression 


ArXi  =  f 


r 


Ki  dx/  =  OaiXQi 


(18) 


SECTION  VI 

THEORETICAL  PREDICTIONS  FOR  THE  SOLAR-INDUCED  ABSORPTION 
PROPERTIES  OF  A  HYPOTHETICAL  THERMAL  CONTROL  MATERIAL 


A  simple  example  will  now  be  given  to  illustrate  the  application  of 
the  mathematical  model  just  described.  Assume  a  hypothetical  mate¬ 
rial  to  be  exposed  to  a  UV  irradiance  whose  magnitude  and  spectral 
distribution  is  equal  to  that  of  the  sun  at  one  astronomical  unit  (Ref.  9). 
During  this  exposure  the  material  is  located  in  a  high  vacuum  and  main¬ 
tained  at  a  constant  temperature  of  approximately  300°K.  The  spectral 
(ultraviolet  and  visible)  absorption  coefficient  of  this  material  prior  to 
exposure  is  assumed  to  be  identical  to  the  one  shown  in  Fig.  4.  Note 
that  €y,  the  absorption  coefficient  in  the  damaging  wavelength  interval 
(0  <  y  <  385  nm)  is  taken  constant  at  a  value  of  3  x  10^  cm-1.  In  turn. 
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R-y,  the  spectral  reflectance  of  the  material  in  this  interval,  is  assumed 
constant  at  a  value  of  0.  07. 


6.1  CONCENTRATION  OF  THE  ith  DEFECT  SPECIE,  N- 

Equation  (9)  is  used  to  calculate  the  concentration  distribution  of  the 
ith  defect  specie.  This  specie  is  assumed  to  have  a  diffusion  coefficient 
in  the  parent  material  of  Di  =  1.  5  x  10" cm^/sec.  The  reference  UV 
spectral  irradiance  in  Eq.  (9),  Ir>  y,  is  taken  to  be  that  of  the  sun,  Is>  y, 
and  the  reference  irradiation  time,  tr,  is  chosen  as  one  second.  For 
these  reference  conditions,  the  spectrally  weighted  (for  solar  irradiance 
Is>  y)  ith  defect  quantum  yield 

U.ydy  (19) 

is  assumed  to  be  10”4  defects  per  photon.  Figure  5  presents  the  defect 
concentration  distribution  in  the  parent  material  with  irradiation  time 
shown  as  a  parameter.  It  is  seen  that  adjacent  to  the  surface  of  the 
material  there  is  a  layer  of  high  and  essentially  uniform  concentration. 
As  a  result  of  diffusion,  the  width  of  this  layer  increases  slowly  with 
time.  Note  that  the  defect  concentration  increases  rapidly  with  irradia¬ 
tion  time  until  a  limiting  value  is  reached. 

Figure  6  presents  the  time  dependence  of  the  defect  concentration 
at  various  depths  in  the  material.  The  curve  labeled  x  =  0  is  the  defect 
surface  concentration,  Ni0,  as  given  by  Eq.  (10).  Note  how  rapidly  the 
surface  concentration  increases  during  the  early  stages  of  irradiation. 
After  this  large  initial  increase,  the  surface  concentration  is  essentially 
constant  with  time.  Similar  behavior  is  observed  for  the  defect  concen¬ 
tration  at  a  depth  of  1.  0  micron.  At  a  depth  of  10  microns,  the  rate  of 
defect  concentration  increase  is  much  less  during  the  initial  stages  of 
radiation.  However,  after  this  initial  radiation  stage,  the  rate  of  con¬ 
centration  increase  at  10  microns  is  greater  than  that  of  the  depths 
closer  to  the  surface  (x  =  0,  1.  0^t).  Hence,  the  defect  concentration  at 
all  three  depths  will  eventually  become  equal. 


kis  - 


r 


kjyr  Is,y  dy 


/f 


6.2  NUMBER  OF  ith  DEFECT  SPECIE  PER  UNIT  IRRADIATED  AREA 

In  Fig.  7,  the  number  of  ith -type  defects,  Q^,  present  per  unit 
irradiated  area  from  the  surface  to  depth  x  is  displayed  as  a  function  of 
depth  into  the  material  with  irradiation  time  as  a  parameter.  This 
figure  shows  that  for  a  given  irradiation  time,  Qi  increases  sharply 
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with  depth  into  the  material  until  a  certain  depth  is  reached.  At  this 
depth,  Qi  becomes  constant  and  remains  so  for  greater  depths.  This 
constant  value  of  represents  the  total  number  of  ith-type  defects 
created  in  the  material  per  unit  irradiated  area  during  the  specified 
irradiation  time.  The  depth  at  which  this  constant  value  of  is 
achieved  increases  with  irradiation  time  because  of  defect  diffusion. 
Note  that  for  a  layer  of  given  depth,  Qj  increases  sharply  with  irradia¬ 
tion  time  until  some  limiting  or  saturation  value  is  attained. 

6.3  INDUCED  ABSORPTION  COEFFICIENT,  A*Ai 

Figure  8  presents  the  increase  in  the  absorption  coefficient  of  the 
parent  material  at  a  depth  of  one  micron  caused  by  the  presence  of 
ith-type  defects.  Irradiation  time  is  shown  as  a  parameter,  and  the 
absorption  coefficient  of  the  material  prior  to  irradiation  is  also  indi¬ 
cated.  In  computing  Ae^,  Eqs.  (9)  and  (15)  were  used.  For  illustra¬ 
tive  purposes,  the  ith-specie  absorption  cross  section,  ^aiA »  was  taken 
as  2  x  10" 17  cm2  and  assumed  to  be  independent  of  wavelength.  From 
Fig.  8,  it  is  seen  that  the  induced  absorption  coefficient,  Ae^,  in¬ 
creases  rapidly  with  irradiation  time  and  at  approximately  10  hr  equals 
the  initial  absorption  coefficient  of  the  parent  material  in  the  visible 
wavelength  range,  2  x  101  cm'*.  At  about  2000  hr,  the  induced  absorp¬ 
tion  coefficient  reaches  a  saturation  value  of  approximately  10  times  the 
visible  absorption  coefficient  of  the  parent  material  prior  to  irradiation. 
However,  in  the  damaging  wavelength  interval  (0,  385  nm),  the  induced 
absorption  coefficient  is  less  than  1/1000  of  the  absorption  coefficient 
of  the  parent  material  prior  to  irradiation.  Thus,  A e-yj/e-y  is  negligible 
and  the  change  in  the  absorption  coefficient  of  the  parent  material  in  the 
damaging  wavelength  range  can  be  ignored.  This,  of  course,  validates 
the  assumption  made  previously  about  e.y  not  varying  with  irradiation 
time. 


6.4  INDUCED  ABSORPTION  OPTICAL  THICKNESS,  A rAj 

Figure  9  displays  Ar^,  the  increase  in  absorption  optical  thickness 
of  the  parent  material  caused  by  the  presence  of  ith-type  defect  specie. 

It  is  presented  as  a  function  of  irradiation  time,  with  depth  into  the 
material  taken  as  a  parameter.  Equations  (12)  and  (18)  were  used  in 
evaluating  At^j.  As  shown  in  Fig.  9,  Ar>^  increases  very  rapidly 
during  the  first  few  hours  of  irradiation.  After  this  period,  At^  con¬ 
tinues  to  increase  with  irradiation  time  but  at  a  much  slower  rate.  Most 
all  of  the  induced  absorption  optical  thickness  At^  occurs  within 
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50  microns  of  the  surface.  For  irradiation  times  of  2000  hr  and  greater, 
this  value  of  At>^  is  approximately  0.  2. 

SECTION  VII 

USE  OF  THE  Nj  RELATION  TO  MATHEMATICALLY  DEFINE  THE  ENVIRONMENTAL 
PARAMETERS  INVOLVED  IN  SOLAR  RADIATION  TESTING  OF  THERMAL  CONTROL  COATINGS 


As  mentioned  in  the  introduction,  one  of  the  objectives  of  this  study 
is  to  mathematically  define  the  space  environment  conditions  which 
should  be  simulated  in  laboratory  tests  of  solar  radiation  effects  on 
thermal  control  coatings.  Since  the  solar-induced  optical  properties 
(Ae\,  At^,  Aa^)*  are  ultimately  a  function  of 

N,  (Q.  -  J**N  jdx) 


Eq.  (9)  will  be  used  to  accomplish  this  objective.  Thus,  with  Ir>  y  taken 
as  the  UV  spectral  irradiance  of  the  sun,  Ig,  y  and  tr  as  some  reference 
solar  irradiation  time,  ts. 


Nj 


(SrM 


yi<T> 

[l  - 


Ry]  kiys(P,T)  [ls>yI 


o  ,yr  f. 


I  ii 

—  1 


'2ey(Ditr 


2(Djl) 


{ 


dy  (20) 


where  k^-yg  is  the  value  of  k^-yr  f°r  the  condition  Ir<  y  =  ISj  y,  tr  =  ts. 


From  Eq.  (20),  it  can  be  seen  that  k^s,  R-y,  and  €y  can  all  depend 
on  the  wavelength  of  the’ damaging  irradiation  and  each  influences  the 
value  of  N^.  Thus,  in  laboratory  tests  of  solar  degradation  effects  on 
thermal  control  coatings,  the  spectral  distribution  of  the  irradiance  on 
the  surface,  I0>  y,  should  closely  simulate  that  of  the  sun  for  the  wave¬ 
length  range  (0,  7i).  This  means  that  the  resultant  spectral  distribution 
of  the  energy  source  and  the  intervening  optics  (e.  g. ,  windows  and 
filters)  between  it  and  the  test  surface  should  closely  approximate  that 
of  the  sun  over  the  interval  (0,  Yi)- 


It  is  also  seen  in  Eq.  (20)  that  the  magnitude  of  the  spectral  irradi¬ 
ance  1Q<  y  influences  the  value  of  N^<  Thus,  for  solar  radiation  testing 
of  thermal  control  coatings,  both  the  magnitude  and  spectral  distribu¬ 
tion  of  the  irradiance  on  the  surface  should  be  nearly  equal  to  that  of  the 


*Aa?^  is  defined  as  the  induced  spectral  absorptance  of  the  parent 
material. 
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sun  in  the  wavelength  interval  (0,  7j).  The  ramifications  of  this  conclu¬ 
sion  for  accelerated  solar  radiation  testing  are  discussed  in  the  next 
section. 

Equation  (20)  also  shows  that  Ni  can  depend  on  the  temperature  of 
the  material  during  irradiation  because  Di,  krys*  anc*  the  wavelength  at 
which  the  material  absorption  edge  occurs,  y^,  all  vary  with  tempera¬ 
ture.  Thus,  the  material  temperature  should  be  controlled  during  a 
solar  radiation  test  and  maintained  at  a  value  near  that  expected  in 
flight . 

Finally,  it  is  seen  that  Nj  of  Eq.  (20)  can  also  be  influenced  by  the 
pressure  of  the  vacuum  environment  in  which  it  is  located.  As  was 
discussed  earlier,  this  is  because  of  the  potential  dependence  of  lq.y  on 
the  concentration  level  of  gases  in  the  surrounding  environment.  Thus, 
solar  radiation  testing  of  thermal  control  coatings  must  be  performed 
in  a  high  vacuum  and  measurements  of  Ae\,  Ar^,  and  must  be 
made  in  place. 


SECTION  VIII 

VALIDITY  OF  THE  TIME-IRRADIANCE  RECIPROCITY  PRINCIPLE  USED  IN 
ACCELERATED  SOLAR  RADIATION  TESTING 


The  discussion  in  the  previous  section  on  solar  radiation  testing  of 
thermal  control  coatings  indicates  that  to  obtain  reliable  data  for  Aa?x, 
both  the  magnitude  and  spectral  distribution  of  the  irradiance  on  the 
surface  should  closely  approximate  that  of  the  sun  in  the  wavelength 
interval  (0,  y^).  This  conclusion,  of  course,  questions  the  reliability  of 
data  obtained  in  accelerated  solar  radiation  testing  where,  to  save  time, 
it  is  assumed  that  the  same  test  results  for  A a\  can  be  obtained  by  sub¬ 
stituting  high -intensity  UV  source  irradiance  for  irradiation  time.  The 
validity  of  this  time-irradiance  reciprocity  principle  is  investigated  in 
the  following  paragraphs  by  use  of  the  equations  developed  earlier  in 
the  study. 

Since  the  induced  spectral  absorptance  of  the  parent  material 
caused  by  the  presence  of  defect  specie  i,  Aa^i,  is  a  function  of  the  in¬ 
duced  absorption  coefficient  Ae^i  =  f(Nj)  and  the  induced  absorption 
optical  thickness  At^  =  f(Q^),  Eqs.  (15),  (9),  (18),  and  (12)  will  be  used 
to  accomplish  the  aforementioned  investigation.  It  will  be  assumed  that 
if  the  reciprocity  principle  does  not  hold  for  Ae^j(Nj)  and  At>^(Q^),  it 
will  not  hold  for  Aarxi-  The  properties  of  the  parent  material  and  the 
ith  defect  specie  will  be  taken  the  same  as  those  used  in  the  previous 
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illustrative  example.  Ultraviolet  irradiance  levels  of  1  and  10  "Suns" 
will  be  used.  The  spectral  distribution  of  both  irradiances  in  the  wave¬ 
length  interval  (0,  7^  =  385  nm)  is  assumed  to  be  that  of  the  sun  (Ref.  9). 
Figure  10  presents  the  surface  (x  =  0)  value  of  the  induced  absorption 
coefficient  Aexi  as  a  function  of  radiant  exposure  in  Equivalent  Sun 
Hours  (ESH)  with  UV  irradiance  level  in  "Suns"  shown  as  a  parameter. 
For  a  given  radiant  exposure,  Ae^  is  seen  to  be  significantly  greater 
at  an  irradiance  of  10  "Suns"  than  at  1  "Sun"  irradiance.  If  the  time- 
irradiance  reciprocity  principle  held,  there  would  be  no  effect  of 
irradiance  level  on  Aexi,  and  only  a  single  curve  equal  to  the  one  "Sun" 
value  of  Ae^  would  appear  in  Fig.  10. 

Figure  11  displays  the  value  of  Aexi  at  a  depth  of  1  micron  in  the 
material.  It  is  seen  that  for  a  radiant  exposure  less  than  60  ESH,  Aexi 
at  1  "Sun"  irradiance  is  greater  than  Aexi  at  a  UV  irradiance  level  of 
10  "Suns.  "  However,  for  a  radiant  exposure  greater  than  60  ESH,  Aexi 
at  10  "Suns"  exceeds  Aexi  at  1  "Sun",  which  was  the  expected  result  for 
all  values  of  ESH.  This  weird  behavior  of  Aexi  with  irradiance  and 
radiant  exposure  in  ESH  is  a  result  of  defect  diffusion.  It  can  be  ex¬ 
plained  by  noting  that  for  an  equivalent  number  of  sun  hours  of  radiant 
exposure  a  defect  generated  by  1  "Sun"  irradiance  can  diffuse  into  the 
material  10  times  as  long  as  a  defect  created  by  an  irradiance  of 
10  "Suns.  "  Thus,  for  depths  in  the  material  which  are  somewhat 
greater  than  the  penetration  depth  of  the  damaging  irradiation,  the 
defect  concentration  Nj  and  hence  Ae^  will  be  larger  for  the  lower  ir¬ 
radiance  level  until  some  particular  ESH  of  radiant  exposure  is  reached. 

In  Fig.  12,  the  induced  absorption  optical  thickness  Ar^  for  a  layer 
of  10  microns  depth  is  presented  as  a  function  of  radiant  exposure  in 
ESH.  The  results  for  UV  irradiance  levels  of  1  and  10  "Suns"  are  pre¬ 
sented.  It  is  seen  that  Atxj  is  larger  at  the  10  "Sun"  irradiance  for 
values  of  ESH  >  200.  For  ESH  <  200,  At^  is  the  same  for  both  irradi¬ 
ance  levels.  The  results  shown  in  this  figure.  Fig.  12,  in  conjunction 
with  those  presented  in  Figs.  10  and  11,  indicate  that  the  time -irradiance 
reciprocity  principle  does  not  hold  for  materials  with  optical  properties 
that  change  according  to  the  photoreaction-diffusion  model  of  solar 
degradation. 


SECTION  IX 
CONCLUSIONS 


A  partial  differential  equation  has  been  derived  for  the  concentra¬ 
tion  of  the  ith  defect  specie  generated  in  thermal  control  materials  by  a 
photoreaction-diffusion  mechanism  of  solar  degradation.  This  equation 
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has  been  solved  to  yield  an  analytical  solution  for  the  defect  specie  con¬ 
centration,  N^.  An  analytical  expression  has  also  been  obtained  for  the 
number  Qj  of  ith-type  defects  present  in  a  layer  of  arbitrary  depth  and 
unit  irradiated  area.  Thus,  a  theoretical  model  for  solar  degradation 
of  thermal  control  materials  can  be  formulated  in  mathematical  terms 
and  analytical  solutions  can  be  obtained  for  both  the  concentration  and 
•number  of  ith  defects  generated  by  solar  UV  radiation.  These  solutions 
may  be  used  to  calculate  both  the  concentration  and  number  of  defects 
created  in  thermal  control  materials  as  a  function  of  irradiation  time 
and  spectral  irradiance.  This  is  a  significant  step  toward  theoretically 
predicting  the  solar-induced  optical  absorptance  change  of  spacecraft 
surfaces  as  a  function  of  time  since  the  induced  absorptance  change  is 
ultimately  a  function  of  defect  concentration. 

The  solutions  for  Nj  and  Qj  can  also  be  used  to  better  define  the 
environmental  parameters  that  must  be  controlled  during  solar  radia¬ 
tion  testing  of  thermal  control  coatings.  For  example,  the  equations 
for  and  show  that  both  the  magnitude  and  spectral  distribution  of 
the  solar  simulator  irradiance  on  the  coating  surface  should  be  nearly 
equal  to  that  of  the  sun  in  the  wavelength  interval  (0,  y^).  In  addition, 
these  equations  indicate  that  the  temperature  of  the  irradiated  material 
should  be  controlled  during  a  solar  radiation  test  and  maintained  at  a 
value  near  that  expected  in  flight.  Also,  the  relations  for  Nf  and 
show  that  solar  radiation  testing  of  thermal  control  coatings  must  be 
performed  in  a  high  vacuum,  and  measurement  of  radiative  properties 
should  be  made  in  place.  Finally,  these  solutions  indicate  that  the 
time -irradiance  reciprocity  principle  used  in  accelerated  solar  testing 
is  invalid  for  materials  which  degrade  according  to  the  photoreaction- 
diffusion  model. 
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x  =  0 


Fig.  1  Mathematical  Model  far  Solar  Degradation  a(  Thermal  Control  Materials 
by  Photoreaction-Diffusion  Damage  Mechanism 


19 


AEDC-TR-68-175 


20 


AEDC-TR-68.175 


21 


Induced  Absorption  Coefficient,  Ae^j,  atx  ■  1  micron  -  cm 
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Fig.  8  Change  in  Parent  Material  Absorption  Coefficient  at  x  =  1.0  Micron 
for  an  Ultraviolet  Irradiance  equal  to  that  of  the  Sun 
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Induced  Absorption  Optical  Thickness,  dimensionless 
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Fig.  9  Change  in  Parent  Material  Absorption  Optical  Thickness  for  an 
Ultraviolet  Irradiance  Equal  to  that  of  the  Sun 
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Fig.  11  Effect  of  Ultraviolet  Irradiance  Level  an  the  Increase  in  Parent  Material 
Absorption  Coefficient  ata  Depth  of  1.0  Micron,  (x  =  lp) 


Induced  Absorption  Optical  Thickness,  Ar^j,  dimensionless 
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Fig.  12  Effect  of  Ultraviolet  Irrodiance  Level  on  the  Increase  in  Parent  Material 
Aborption  Optical  Thickness  at  a  Depth  of  10  Microns,  Ar\i  (x  =  10/i) 
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